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A useful and universal formula for the expectation value of the radial operator in the presence of 
the Aharonov-Bohm flux and the Coulomb Field is established. We find that the expectation value 
(r x ) (-00 < A < oo) is greatly affected due to the non-local effect of the magnetic flux although the 
Aharonov-Bohm flux does not have any dynamical significance in classical mechanics. In particular, 
the quantum fluctuation increases in the presence of the magnetic flux due to the Aharonov-Bohm 
effect. In addition, the Virial theory in quantum mechanics is also constructed for the spherically 
symmetric system under the Aharonov-Bohm effect. 

PACS: 03.65.Sq 



I. INTRODUCTION 

The calculation of the expectation values (r*) with different power A of the radial operator f is important in the 
Hellmann-Feynman theory [Q,^), the interactions of the molecular theory ||, and the quantum chemistry j^J^,^,^|. 
On the other hand, the Aharonov-Bohm (AB) effect, a topological non-local physical significance at the quantum 
level, has shed light on the understanding of the phenomenon of the fractional quantum Hall effect [pp|, ^0| , pd[ | , the 
superconductivity [[y]jlj], and the repulsive Bose gases [Ol in the last 20 years. It is known that the global influence 
of the AB effect affects all charged particle's systems [p~4| . In this paper, we will derive a recurrence formula for 
the expectation value (f*) (— oo < A < oo) of the radial operator f in the presence of the Coulomb field under 
the Aharonov-Bohm effect (ABC). Moreover, the Virial theory in quantum mechanics will also be generalized to the 
system with spherical symmetry under the AB flux. The general effect of the AB flux to the expectation value will 
also be discussed in this paper. 

The fixed-energy bare Green's function G°(r, r' ; E) for a charged particle with mass m propagating from r' to r 
satisfies the Schrodinger equation 



E-H Q (r,-V) 

i 



G°(r, y'-E) 



(1.1) 



Here the Hamiltonian of the system is given by Hq — — h 2 V 2 /2m + V(r). The angular decomposition of the Green's 
function can be written as 



G°(r,r';E) =£ £ G°(r/; E)Y lk (9, ^1(6', </), 



(1.2) 



2=0 k=-l 



i n th e spherically symmetric system with the well-known spherical harmonics. Hence the left hand side of the Eq. 
(1.1) can be reduced to the following form: 



oo I 



1=0 k=~l 



2 d \ + i)h 2 



2m V dr 2 r dr 



2mr 2 



V(r) 



xG?(r,r';£)H t (^)yi;(«>'). 



(1.3) 
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II. RADIAL WAVE EQUATION FOR A CHARGED PARTICLE IN A MAGNETIC FIELD 



For a charged particle in a magnetic field, the charged Green's function G is related to the bare Green's function 
G° by the following equation: 



G(r, r'; E) = G°(r, r'; E)e^ £ A ^ d \ 



(2.1) 



with a globally path-dependent nonintegrable phase factor |15|,[l6| given above. Here the vector potential A(r) is used 
to represent the contribution from the magnetic field. For the Aharonov-Bohm magnetic flux under consideration, 
the vector potential can be written as 



A(x) 



2" p e ¥> 27rp e V 



(P<e) 
(P>e) 



(2.2) 



Here the two-dimensional radial length square is defined as p 2 



x 2 + y 2 . Moreover, e v is the unit vector of the 
coordinate (p and e is the radius of region where the magnetic field exists. Hence the total magnetic flux is given by $ = 
ire 2 B. Note that the associated magnetic field lines are confined inside a tube, with radius e, along the z-axis. Along 
the region free of the magnetic field, the path-dependent nonintegrable phase factor is given by exp{— ipo Jp dr'ipij 1 )}. 
Here we have used the subscript P to represent the path-dependent nature of the phase factor. In addition, we have 
denoted ^(r') = dip/dr'. Moreover, p,o = —2eg/hc is a dimensionless numerical factor defined by $ — Aitg. The 
minus sign we adopted is a matter of convention. According to the discussion in Ref. Jl6[ , only phase factors with 
closed-loop contour are considered where the description of electromagnetic phenomenon are complete. Hence, we 
have 



(2.3) 



with integer values n corresponding to the winding number. The magnetic interaction is therefore a purely topological 
phenomenon. Therefore the nonintegrable phase factor becomes exp{— i/i (2ri7r)}. With the help of equality between 

the associated Legendre polynomial Pj?{z) and the Jacobi function P„ (z) 0,^8), we find that 



P, fc (cos< 



(-1) 



fe r(z 



fc + i) 



cos -sin-) p/*f 



(cos (9). 



r(z + i) 

Therefore the angular part of the Green's function in the expression (|l.3|) can be shown to be 



(2.4) 



£ y ft (My,S(*V) = £ 2 4 + 1 r a + 1 1 1! p " (cos g)p " (cos *'K fe(y - y,) 



E 

k=-l 



2l + lT(l-k + l)T(l + k+l) 



4tt 



r 2 (z + i) 



cos — cos — sin — sin — 
2 2 2 2 



P^'P (cos e)P, ( *^ (cos e')e lk ^') . 



(2.5) 



To include the nonintegrable phase fac tor due to the AB effect, we can rename the index I into q related by the 
relation I — k = q. As a result, Eq. ( |l.3| ) can be written as 



E - 



OO OO 

E E 

q—0 k— — oC' 





( d 2 


2 d \ 


{q + k)(q + k + l)h 2 


2m 


\dr 2 


r dr J 


2mr 2 



V(r) 



xG° q+k (r/;E) 



2(q + k) + IT (q + 1)T {q + 2k + 1) 



4tt 



T 2 (q + k + 1) 



cos - cos — sin - sm — 
2 2 2 2 



2 



x p( fc ' fc ) (cos 8)Pj l k ' k '> (cos e')e ik ^-'p') . 



(2.6) 



In addition, the nonintegrable phase factor exp{— (2mr)} can thus be included with the help of the Poisson's 
summation formula (p. 124, Jl9||) 



£ /(*) = / dyJ2 



a 2nnyi 



k= — oo 



Hence, the expression (2.6) can be written as 



v q=0 k—~oc 



2m V dr 2 r dr 



(q + z)(q + z + l)fi 2 



2mr 2 



V(r) 



(2.7) 



xG q+z (r/;E) 



2(q 



ir(g + l)r(g + 2z + l) 



47T 



' \ z 



T 2 (q + z + l) 
x p( z < z > (cos 8)Pj 1 z ^ (cos ^) e *(»-rt))(v+2far-(p') . 



cos — cos — sin — sin — 
2 2 2 2 



(2.8) 



Here the superscript in G° q+k has been suppressed to reflect the inclusion of the AB effect. The summation over all 



indices k forces z = pbQ modulo an arbitrary integral number. Therefore, one has 

(q+ \k + n \)(q+ |fc + Mo| + l)^ 1 



OO OO 



E -E E 

^ q—0 k— — oo 



[-— 


f d 2 


2 d \ 


2m 


\dr 2 


r dr J 



2mr 2 



V(r) 



vr , , F , j [2( g +|fc + M o|) + l]r(g + l)r(2|fc + ^o|+ g + l) \ ifcf /) 



x (cos 6»/2 cos 6»'/2 sin 6»/2 sin 6»72) |fe+Ml P (Z (|A;+ ^ Mfc+ ^ l) (cos6l)P ( 



(|fc+/i |,|fc+A»ol) 



{cos 9'). 



(2.9) 



Note that the effect of the AB flux to the radial Green's function is to replace the integer quantum number I with 
the fractional quantum number q + | k + //q | . Analogously the same procedure can be applied to the delta function 
<5 3 (r — r') in the rhs of the Eq. (1.1) with the help of the following solid angle representation of the 5 function: 



OO I 



*(fi-fi')=E E Y lk {0^)YC k {0',v')- 



(2.10) 



1=0 k=-l 

Therefore, for the set of the fixed quantum numbers (q, k) one can show that the radial Green's function satisfies 



E 



_ I A. 

2m V dr 2 r dr 



{q+\k + ^\){q+\k + n \ + l)Tl 2 



2mr 2 



V(r) 



xG q+lk+fiol (r,r';E)=S(r-r') 
Hence, the corresponding radial wave equation reads 

ft 2 a{a + 1) 



h 2 d 2 

Ira dr z 



E - V(r) + 



u a (r) = 0, 



(2.11) 



(2.12) 



2m r 2 

where we have set a = q + \k + [1q\, and u a (r) — rR na {r). It is clear that R na satisfies the spherical Bessel equation 



d 2 , 2 d | ( ? a(a+l) 
dr 2 r dr \ r 2 



Rna(r)=0 



(2.13) 
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with the definition k = \J 2mE /ft 2 and the reduced potential U(r) = 2mV (r) / ft 2 . For simplicity, we have written 

Rna(r) instead of R n ,q,k{ r ) m which each set (n,q, k) denote a quantum state. Hence the AB effect reflects itself by 
the coupling to the angular momentum in the radial Green's function, which turns the integer quantum number to a 
fractional one. 

In order to derive the expectation value (V A ) in an arbitrary quantum state (n, g, k) of the ABC system, one notes 
that the attractive potential of the Coulomb field is given by V(r) = —Ze 2 /r in which Ze is the total charge at the 
center of source and — e representing the charge of the electron. One can show that the exact solution of the energy 
spectra in this case is given by 20 23| 



E 



n : q : k 



2clq [n 



where do = ft 2 /me 2 is the Bohr ra dius, and the ranges of the quantum number are n, q 
With the help of the result, Eq. fl2.12| ) can be brought to the following form: 



0,1,2, 



(2.14) 



and — oo < k < oo. 



dr 2 



u(r) 



2Z a(a + l) / Z 
aor r 2 \nao 



u{r) = 0. 



(2.15) 



III. RECURRENCE FORMULA FOR THE EXPECTATION VALUE OF THE RADIAL OPERATOR OF 

THE AHARONOV-BOHM-COULOMB SYSTEM 



For simplicity, we have written u(r) instead of u a (r), and n = [n + q + \ k + /xo| + 1]- For our purpose of calculating 
the diagonal matrix element, we will multiply r x u to both sides of the above equation. Integrating this equation with 
respect to r, J Q ■ ■ ■ dr, one derives 



r x u 



*(a + 1) (r x - 2 ) + — (r^ 1 ) - ( JL V ( r A\ = . 
a \na J 



(3.1) 



With the help of integration by part, the first term yields 



r x u [ ~J^ U ! =: " 



dr 



r x —u + Xr 
dr 



a-i„ 



dr 



u dr 



r"u [ — u I — — r u 
dr 



A(A-l) (r A_ 2 x 



\ I d 
r x 



dr 



u dr. 



One can show that the asymptotic behavior of the wave function is given by |T^ , ^2]| 



r — > 0, u ~ r a+1 

r > OO, U ~ r n e -Zr/na 



(3.2) 



(3.3) 



directly from the asymptotic wave equati on. Therefore, one has the necessary and sufficient condition to guarantee 
the vanishing of the first part of the Eq. (3.2) 



/■ ;/ I — u 



0, r^V 1 



0, 



as long as A > — (2a +1). Therefore 
" A ( A -!) 

— - a(a + 1) 



<r*-»> + ^<r*-i>-(A) <^> = 



—u ] ilr. 

dr 



(3.4) 



(3.5) 



Similarly, let's multiply 2r x+1 du/dr to the both sides of the Eq. (2.15) and then integrate over J °° ■■ - dr. One can 
thus derive the following equations: 
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Jo 



+ dr, 



2r x+1 u 



(3.6) 
(3.7) 



One can also der ive s imilar equations with the power factor A + 1 repl aced by A or A — 1 in above equation. According 
to the equation (^4), the first terms in the rhs of the Eq.s (3 1 6)-(3/7) all vanishes. Hence, we obtain 



(A - l)a(a + 1) (r^) 2A^ (r^) + (A + 1) (A ) <r>) 

a \nao,' 



= (A + 1) / r 
lo 



—u I dr. 

dr 



(3.8) 



Eliminating the term on the rhs of the above equation with the help of the equation (3J5), we finally obtain the general 
recurrence formula for the ABC system 



A 



±1 (r*) - (2A + 1)| (r^) + \ [ { 2a + i f - A 2 J (r^) 



= 0. 



(3.9) 



This formula provides a very convenient tool to calculate the expectation value of the radial operator f x for arbitrary 
power A in any arbitrary quantum state (n, q, k) for the ABC system without facing the complication computing 
the wave functions directly. For example, let's calculate a few leadin g ter ms of the expectation value of the radial 
operator. First of all, one can set A = where the recurrence formula (3.9) gives 



1\ Z Z 

r /n, q ,k " 2fl o (n + q+\k + Mo| + l) 2 ao' 



(3.10) 



Note that the effect of the magnetic flux leads to the decrease of the expectation value of the potential. Note that 
this effect has no correspondence in classical system. Moreover, if we choose A = 1,2, one will instead derive 



z ■ 



^) n , q , k = \[^ 2 -a(a + l)] 
<^W = ¥[l + 5n 2 -3a(a + l)] (f) 



(3.H) 



These results show that the magnetic flux effect will increase the expectation of the radial operator such as the case 
where the quantum states given by (n, q, k) — (0, 0, 0). Indeed, one can show that 



( r )o,o,o ~ 
V /o,o,o 



2|/xof 
wl+i) 2 



' 5 |/iol 



"0 

z ■ 



2|Mor + 7|Mo|+6 



(ff 



(3.12) 



Note that these results reduce to the well-known pure Coulomb system, where = 0, (r) = 3ao/2Z, and 
(r 2 ) Q Q = 3(ao/Z) 2 where the nonlocal effect of the magnetic flux in the physical average becomes manifest in these 
formulae. 



According to the definition leading to the Eq. ( 2.15| ), we have a 



Mo I 



1 corresponding to 



the nodes of the wave function at n 
reduces to Q 



for the modified circular Bohr orbit. In this case, the radial wave function 



Uo, q ,k(r) = rRo t(ltk (r) = r 



(3.13) 



Hence we can calculate the most possible position of the radius r most , defined as the position where the probability 
function |m| 2 attains its maximal, given by <iuo i9i fc(r = r most )/dr = 0. The result indicates that 



h 2 ao 
Z 



(3.14) 



for the most possible position of the radius. It is clear that the presence of the magnetic flux increases the value of 
the most possible position. On the other hand, Eq. (3.11) gives the following radial expectation value, when a is 
equal to a max , 



5 



( r >o,s-i 



r,2 



oo 
Z ' 



(3.15) 



Therefore quantum mechanical average (the expectation value) for the radius of the modified Bohr orbit is bigger 
then the most possible position of the radius. In addition, the position fluctuation is given by 



n 3 n 2 



Z 



(3.16) 



which implies the increase of the quantum fluctuation due to the increase of the magnetic flux. Consequently, one 
can derive the ratio between the fluctuation and the corresponding average 



1 



<r>o, fi -i V^n + T 



(3.17) 



which indicates that the quantum mechanical result tends to approach the Bohr picture of orbit quantization as n 
increases or equivalently the magnetic flux increases. 

If we choose A = — 1, —2, the recurrence relation indicates that 



/ r -3V = ^O/Z , 2 . 

\ /n, q ,k a(a + l){a a /Zy K ' 



1 



n,q,k n 3 a(a + 1/2) (a + 1) \a Q 



(3.18) 



, _ 4 > 3ri 2 — a(a + 1) ( Z 

^ < n 'i' k ~ 2n 5 {a - l/2)(a + l/2)(a + l)(a + 3/2) 



Here we have used the result 



\ /n.q.k n 3 (a + l/2) \a a J 



(3.19) 



(3.20) 



which is the only moment can not be derived from the recurrence Eq. ( |3.9| ). It can, however, be derived from the 
well-known Hellmann-Feynman formula 



8E; 



da 



n,q,k 



dH 



Da 



n,q,k J 7 



where the Hamiltonian operator is given by 



H 



h 2 d 2 h z a(a + 1) Ze 2 
2m dr 2 2m r 2 r 



It hence follows that 



q . k 



Oa 



(3.21) 



(3.22) 



(3.23) 



With the help of the Eq. (2.14), one derives the result shown in Eq. (3.20). This result gives us the ratio information 
of the modified centrifugal potential in quantum mechanics, namely, 



a(a + l)h 
2mr 2 



a(a + l) 
n, q ,k ~ ~(a + l/2)n 



Ef,. 



(3.24) 



Since the total kinetic energy of the system equals < T > n ,q,k= —Efi (see appendix for details), the ratio of the 
centrifugal potential to the total kinetic energy, r c =< V c > / < T >, is 

r c = a(a + 1)/ [(a + 1 /2)n] . 

This result indicates that, for fixed h, the increase of the intensity of the magnetic flux induces the increase of the 
ratio r c and thus the weight of the centrifugal potential. When a is equal to a max = h — 1 corresponding to the 
modified Bohr orbit case, the ratio r c become 



G 



On the other hand, the ratio of the average of the radial kinetic energy to the total kinetic energy, r r =< p 2 /2m > 
j < T >, is given by 



2m 



1 



0,q,k 



2ft- 1 



(3.26) 



at the limit a = a max . Therefore when h 3> 1, the radial kinetic energy will become very small reproducing the 
classical result. 

In addition, one can also discuss the effect of the magnetic flux in the well-known formulae introduced by J. 
Schwingcr plj . It was shown that 



(3.27) 



These results are very useful and of considerably general interest for the study of the bound states in a system with 
a central potential. The first identity relates the s-wave wave function at the origin to the gradient of the potential. 
This identity comes explicitly with the Planck constant ft indicating the significance of the quantum effect. Therefore 
there is no classical correspondence. The second identity is, however, a generalization of the quantum effect of an 
important classical theory 



dV 
dr 



l 2 



f 



(3.28) 



where (/) av represents the average of the periodic function f(t) over a classical period. Eq. ( |3.28] ) states that the 
average of a external force equals to the average of the centripetal force ov er th e classical period. 
In the presence of the magnetic flux, the Hamilton operator is given by (2.13) 



H = - 



2m \ dr' 



2 8_ 

r dr 



a(a + l)h 2 
2mr 2 



V(r). 



Computing the commutator of the operators d/dr and H, one has 



d_ 

dr ' 



H 



d a(a + l)h 2 



r 2 



dV_ 

dr 



(3.29) 



(3.30) 



Evaluating the expectation value with respect to an arbitrary bound state which makes the contribution from the lhs 
vanished, we would have 



dV 
dr 



a{a + l)Ti 2 



1 



m 



LJL\ 

r 2 dr I 



1 d 



r 2 drdn^*^^, 
ml r z dr 



(3.31) 



with dtt the solid angle. Since the mean value in the lhs of the above equation is a real number, so does the rhs of 
the above equation. Hence one has 



dr^*-^-ty 

dr 



9 » 
dr 



d 

dr— 

dr 



1 



(3.32) 



In deriving the above equations, we have used the following facts: (1) the wave function of the bound state vanishes 
at the region r — * oo; (2) in the presence of the magnetic flux, the vanishing of the wave function ^(0) = holds 
at the origin p2| . To prove the second proposition, one notes that the asymptotic form of the wave function can be 
shown to be |E2f 
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a = q + \k + /i | > 0, 



(3.33) 



near the origin in the presence of the magnetic flux. Because that u a (r) — r R a (r ), one has immediately that 
u' a (r) = R a (r) + rR' a (r). Therefore we obtain u' a (0) = R a (0). In addition, Eq. ( |3.33| ) says that <(r) = (a + l)r a 
implying that u' a (0) — 0. This proves the second proposition that ^(0) = 0. 



Accordingly, in a central force system with a magnetic flux, one can generalize the Schwinger theorem (3.27) to the 
following form: 



dV\ 

dr /n,g,k 



a(a + l)h 2 



1 



(3.34) 



n,q,k 



In addition, the s-wave no long er ex ists in the presence of the magnetic flux, there is no corresponding generalization 
of the first identity in the Eq. ( 3.27 ). On the oth er ha nd, we find that the topological effect of the magnetic flux has 
a generalized quantum mechanical shown in Eq. (3.27) although such effect does not have any dynamical significance 
in classical mechanics. 



IV. CONCLUSION 



We have derived a general recurrence formula of the expectation value of the radial operator with different power for 
an Aharonov-Bohm-Coulomb system providing us a convenient method to compute the average of the radial operator 
in any power. With this useful recurrence formula, one can avoid tedious and complicate integration with the wave 
functions. The useful formulae and techniques shown in this paper have many interesting applications in the study 
of many physical systems in the presence of a magnetic field. For example, we have shown the generalized theorem 
of Schwinger for bound state system under the influence of a magnetic field. 

In addition, the two dimensional (2D) system is important in the fractional quantum Hall effect and the high-T c 
superconductivity, we will also generalize our result to the 2D system at the end of this paper. Two dimensional 
central force quantum system is governed by the radial Schrodinger equation |25|] 



2m 



d?_ 
dp 2 



1 d 
pdp 



,,2 



V{p)\R nk = ER nk , 



(4.1) 



where p is the corresponding radial length in 2D, i.e. p — \J 1 x 2 + y 2 . The corresponding total wave functions is given 
by ^(p, ip) = R n k(p)e lkip with the range of the quantum number specified by k — 0, ±1,±2, In addition, the 
energy spectra E depends on |fc| and the wave function R n \ k \ corresponding to the energy level E n \u satisfies the 
orthonormal condition 



dppR* n \k\(p) R n'\k\(p) = $n 



If we perform a transformation R n \k\{p) = Xn\k\{p) I ' \J~P-, Eq- (4.1) will reduce to 



2m 



(fc-l/2)(fe + l/2) 



dp 2 



+ V(P) >Xn\k\(p) = EXn\k\(p)- 



The AB effect can be introduced by the similar procedure leading to the Eqs. (2.7)^( |2~13| ). The result is 



2 r A2 



2m 



(a- 1/2) (a + 1/2) 



dp 2 



+ V(p) \ Xn\k\(p) = EXn\k\(.P), 



(4.2) 



(4.3) 



(4.4) 



where we have defined a = \k + /xo | - Comparing this equation with Eq. ( 2.12 ), one finds that results derived for the 
3D system can be generalized to the 2D system with the parameters a and a related by the following relation 



or simply 



a(a + 1) < — ► (5- l/2)(fi + 1/2), 
a< > (5-1/2). 



(4.5) 
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For example, we find that the energy spectra of the 2D ABC system is given by 

Z 2 e 2 Z 2 e 2 



2a [n+ |fc + /io| + 1/2] 2 a n 2 



(4.6) 



This exact formula agrees with the previous result obtained in Ref. |2^]. In addition, one can also derive some useful 
and important expectation values for the 2D ABC system: 



Z 

P/n.k n 2 



(4.7) 



1 ( Z 



P 2 I ' n,k ^2" \ao 



(4.8) 



P° 



n , k 



z 

n\a{a- l/2)(a+ 1/2) Va^ 



(4.9) 



Generalization to many other cases are straightforward following similar procedures. 
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V. APPENDIX 



In the appendix, we shall prove the Virial theorem for the spherically symmetric system under an AB magne tic 
flux. The Hamiltonian for the spherically symmetric system under an AB magnetic flux is given by the Eq. ( 2.13 ) 



2m Zmr^ 



where the Hermite operator is defined as 



h ( d 1 



(5.1) 



(5.2) 



To obtain the Virial theorem, one notes that the Heisenberg equations of motion for the position and the momentum 
operators are 



dr _ 1 

dt ih 



r.H 



Y7 £t - P r 

= VpH = e Pr — , 
m 



(5.3) 



dp _ 1 

dt ih 



= V r H = 



a(a + l)h 2 dV{r) 
mr 3 dr 



(5.4) 



where e Pr is the unit vector of p r in the momentum space. Therefore, we have 

d 



dt 



(r-p) = 



1 

ih 



rp,H 



dr 



dt P + r dt 



dp _ f>l a(a + l)h 2 r dV(r) 



dr 



(5.5) 



Average with respect to the wave functions, we obtain 







pi a(a + lW\ 1 / dV(r) 

— r 



2m 2mr 2 



dr 



(5.6) 



where we have used the identity 



rp,H 



T = 



0. Let us define the total kinetic energy operator as 
p 2 r a(a + l)h 2 



(5.7) 



2m 2mr 2 

Hence one can derive the following Virial theorem for the spherically symmetric system under an AB magnetic flux 

dV(r 



2(T 



dr 



(5.8) 



In particular, if the potential V(r) is a homogeneous function of degree v, namely V(br) = b u V{r), the Eq. ( |5.8| ) can 
be shown to be 



2(f)=v{V{r)) 



Together with the total energy relation, 



(5.9) 



(5.10) 



it is easily to show that the relations between the averages of the kinetic and the potential energy, and hence the total 
energy 



(5.11) 



' n,q,k 



For an ABC system (y = —1), one has < T > n , q ,k= — -En,q,fc, and < V > n , q ,k= 2E n ^k- 
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